Equilibrium sets in quintom cosmologies: the past asymptotic dynamics 
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In the previous paper [3| was investigated the phase space of quintom cosmologies for a class of 
exponential potentials. This study suggests that the past asymptotic dynamics of such a model can 
be approximated by the dynamics near a hyperbola of critical points. In this paper we obtain a 
normal form expansion near a fixed point located on this equilibrium set. We get computationally 
treatable system up to fourth order. From the structure of the unstable manifold (up to fourth 
order) we see that that the past asymptotic behavior of this model is given by a massless scalar field 
cosmology for an open set of orbits (matching the numerical results given in [3). We complement 
the results discussed there by including the analysis at infinity. Although there exists unbounded 
orbits towards the past, by examining the orbits at infinity, we get that the sources satisfy the 



evolution rates (j> /V 



2m-' 



%ln, with Hj^ approaching zero. 



PACS numbers: 47.10.Fg, 05.45.-a, 02.30.Mv, 98.80.-k, 95.36.+X, 98.80.Jk 



I. INTRODUCTION 

Since 1998, the observational evidence coming from 
type I-a Supernovae (SNIa), Large Scale Structure 
(LSS) formation, and the Cosmic Microwave Background 
(CMBR) Radiation ^ 

for an expanding Universe is more 
and more convincing. The acceleration of the expan- 
sion seems to be fuelled by a yet unknown kind of mat- 
ter named dark energy (see the reviews 0). The most 
promising particle candidates for dark energy are spin 
zero bosons, which in the mathematical formalism of 
Quantum Field Theory are represented by scalar fields. 
Many models of dark energy use a single scalar field and 
two very popular options are the so called phantom field 
and the quintessence field. 

Beyond the theories with a single scalar field, models 
with two fields (quintessence and phantorn) have settled 
out explicitly and named quintom models d, 0, H, 0, H, 
H) El, [3 • The quintom paradigm is a hybrid con- 
struction of a quintessence component, usually modelled 
by a real scalar field that is minimally coupled to grav- 
ity, and a phantom field: a real scalar field -minimally 
coupled to gravity- with negative kinetic energy. Let us 
define the equation of state parameter of any cosmolog- 
ical fluid as w = pressure/density. The simplest model 
of dark energy (vacuum energy or cosmological constant) 
is assumed to have w — —\. A key feature of quintom- 
like behavior is the crossing of the so called phantom 
divide, in which the equation of state parameter crosses 
through the value w = —\. Quintom behavior (i.e., the 
w — —\ crossing) has been investigated in the context 
of h-essence cosmolo gies [1, in the context of holo- 
graphic dark energy [14| ; inspired by string theory Q ; 
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derived from spinor matter [T3|; for arbitrary potentials 
[3 [3 ; using isomorphic models consisting of three cou- 
pled oscillators, one of which carries negative kinetic en- 
ergy (particularly for investigating the dynamical behav- 
ior of massless quintom) jl^. The crossing of the phan- 
tom divide is also possible in the context of scalar tensor 
theories [20l.[2ll . [22 . l23l [2l | as well as in modified theories 
of gravity [25l |. 

The cosmological evolution of quintom model with ex- 
ponential potential has been examined, from the dynami- 
cal systems viewpoint, in [6*1 and The difference be- 
tween [6i] and 0, @] is that in the second case the potential 
considers the interaction between the conventional scalar 
field and the phantom field. In the case in which the 
interaction term dominates against the mixed terms of 
the potential, was studied. It was proven there, as a dif- 
ference with the results in , that the existence of scaling 
attractors (in which the energy density of the quintom 
field and the energy density of DM are proportional) ex- 
cludes the existence of phantom attractors. Some of this 
results were extended in p^ . for arbitrary potentials. 

In 7] were investigated curves of (non-hyperbolic) crit- 
ical points C±, and other types of hyperbolic critical 
points from the dynamical systems perspective. We will 
focus here on the curves C±. These curves are associated 
with the past dynamics of quintom models [7] ■ This fact 
was supported by numerical integrations. Although can 
be argue that these curves are non-physical because the 
value w = \, we believe that the investigation of the past 
dynamics of such system is important for completeness. 
This involve the analysis of sets of non-hyperbolic critical 
points. 

When a vector field x' = X(x) in M", admits non- 
isolated critical points, for instance a curve C of criti- 
cal points, then the matrix of derivatives DX(x), when 
evaluated on each point in the equilibrium set C, has 
necessarily one zero eigenvalue. Hence, each point in C 
is non- hyperbolic. Then, linealization technique fails to 
be applied. Although the equilibrium points in C are 
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non- hyperbolic, one is able to apply the Invariant Man- 
ifold Theorem [2^, [S^, [1^ to guaranteed conditions for 
the existence of stable and unstable manifolds. Suppose 
that each point in the equilibrium set C, assumed to be 
a curve, has an unstable manifold of dimension n„, the 
union of these manifolds forms an (ri„ + l)-dimensional 
set whose orbits approach a point of C as t —oo. If 
Uu = n — 1 we say that C is a local source [1^ . If the 
equilibrium set has only one zero eigenvalue and the oth- 
ers with non-zero real parts, then the equilibrium set is 
called normally hyperbolic [HI. In such a case one is 
able to obtain useful information about the stability of 
the critical set by examining the signs of the other, non- 
null, eigenvalues, as we will see later. However, if one is 
interested in the construction of center manifolds, or at 
least approximated ones, if conditions for their existence 
are fulfilled (for references on center manifolds see [27j . 
chapter 4; [s^], [Is'Kb), chapter 18), we require to use 
other tools, e.g., normal forms. 

For the construction of (at least, approximated) in- 
variant manifolds for critical points, we can transform a 
dynamical system defined in the neighborhood of a criti- 
cal point to a non-linear dynamical system in its simpler 
form. The theory of Normal Forms (NF) can be applied 
in order to do so, and in most occasions, terms of order 
higher than two in the Taylor expansion are required. 
For the construction of normal forms for vector fields in 
M" we will follow the approach in [l^l, chapter 2. NF 
calculations consist of two stages: first, the construction 
of the normal forms in which case the nonlinear terms 
take their simplest form and second, the determination of 
the topological types of the fixed points from the normal 
forms. Then by using techniques described in 28] (b), we 
will able to construct approximated invariant manifolds 
up to the desired order in the vector norm. 

In this paper we try and investigate Normal Forms 
expansions of the system in [3] near an hyperbolae of 
critical points. Our main purpose is to devise approx- 
imated unstable and center manifolds. We develop the 
normal expansion (with undetermined coefficients) up to 
an arbitrary order. The system resulting by truncating 
error terms admits a solution given in quadratures, which 
in general cannot be solved explicitly. By integrating the 
resulting system, one is able to construct (theoretical) ap- 
proximated unstable and center manifolds. The problem 
to determine all the coefficient values, up to a prescribed 
order, if very difficult. Our computing power allow us to 
determine the coefficients values up to fourth order. 

The paper has been organized as follows: in section ITTl 
we make a general description of the model under study. 
In section IIIII we recall and improve some results from 
reference jTj. To make the paper self-contained we ex- 
plicitly state in section [A] the main techniques of the NF 
expansions by following the approach in [13] (see also 
[2^(b), chapter 19). We apply these techniques to the 
quintom scalar field solution in section IIVI (which we ar- 
gue are associated to the past asymptotic dynamics). In 
section |V] we offer the analysis at infinity obtaining suffi- 



cient conditions for the existence of local sources at this 
regime. In section IVTl wc offer our main conclusions. 



II. THE QUINTOM MODEL 

In the following discussion on the quintom phase space 
analysis we restrict ourselves to the two-field quintom 
model, with a Lagrangian: 



2 ^ 



-d^cpd^ip-V{^,ip), (1) 



and we include, also, ordinary matter (a comoving perfect 
fluid) in the gravitational action. As in [3] we consider 
here the cfective two-fleld potential 



(2) 



where the scalar fleld (j) represents quintessence and ip 
represents a phantom field. For simplicity, we assume 
m > and n > 0. 

The geometry is given by the Friedmann-Robertson- 
Walker (FRW) metric (in spherical coordinates): 

ds^ = ^dt^ + a{t) (dr^ + r^de^ + sin^ OM^) . (3) 

The field equations derived from ([3]), are 



H'^-lU^-^^)-lV-yM = 0, (4) 



(02_^2^ + iy_ip^, (5) 

Pm = -3-ffpM, (6) 
^ + 3H4. - V6mV = 0, (7) 
(f + 3Hip + VGuV = 0, (8) 

where H — denotes de Hubble expansion scalar. 
The dot denotes derivative with respect the time t. 



III. SOME RESULTS FOR THE FLAT FRW 
CASE 

In [3] it was studied an homogenous and isotropic uni- 
verse, having dust (DM) and quintom DE using the stan- 
dards tools of the theory of dynamical systems. There are 
introduced the normalized variables: {x^, x^, y), defined 

by 



V 



(9) 



They are related through the Friedman equation (U by 



Pm <- 1 
3H^ ^■ 



The dynamics in a phase space is governed by the vec- 
tor field (or differential equation) [3]: 
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x'^^- (3my2 + _ 2)x^) (10) 
< = -J(3V-(g-2)x^) (11) 
y' = ^{l + q-'i{mx^+nx^))y (12) 

defined in the phase space given by 

* = {x = (x^, x^, y):Q<xl~xl + y^ < 1}. (13) 

Here the prime denotes differentiation with respect to 
a new time variable r — log , where a is the scale factor 
of the space-time. The deceleration factor q = —aa/a?' 
can be written as 

g = i (3 - 4 - y^) + 1) . (14) 

A. Linear analysis of non-isolated critical points 

In this section we want to investigate, in more de- 
tail, the behavior of the two curves (hyperbolae) of 
non-isolated fixed points C± with coordinates — 

±^l + a:*^, x^ = x* , y = 0. 

From the physical viewpoint, the curves C± represent 
solutions in which the contribution of matter and the 
potential energy to the total energy density is negligible. 
These solutions are therefore of stiff- fluid type {w — 1), 
which in turn correspond to a decelerating universe. On 
the other hand, the curves C± are local sources, and 
therefore it is unlikely that they can represent the final 
stage in the evolution of our Universe. Thus, by analyz- 
ing the sign of the real part of the normally-hyperbolic 
curves C± we get the following results (we are assuming 
TO > and 71 > 0): 

1. If TO < n, C-i- contains an infinite arc parameter- 
ized by Xl such that xt < -n-mVl-rn-^+ v^ ^^g^^ j 

a local source. C_ contains an infinite arc parame- 
terized by such that x*^ < that 
is a local source. 

2. If TO = n, C+ contains an infinite arc parameterized 

1 _ ^ 

by X* such that x* < that is a local source. 

All of C_ is a local source. 

3. If TO > n there are two possibilities 

(a) If TO^ — n'^ < 1, all of C_ is a local source. 
A finite arc of C+ parameterized by x* such 

is a local source. 

(b) If — r? > 1, no part of C+ is a local source 
and all of C_ is a local source. 



In conclusion, we have that C_ is always a local source 
provided m > n > 0. But, we said nothing about its 
global stability. To answer this question, one needs to 
apply more sophisticated tools such as the calculations 
of Normal Forms (we submit the reader to section HVl for 
such an investigation). 

B. Heteroclinic sequences 

From the local stability properties of C±, O, P and T 
and supported by several numerical integrations Q it 
was possible to identify heteroclinic sequences (see table 
Ufor additional information concerning the critical points 
0,C±,r,P reported in [3]). 

• Case i a) If < TO < TO, the point P is an stable 
node and T does not exists. For this conditions 
there exist an heteroclinic sequence of type /C_ — > 
O — > P or of type IC+ — > O — > P, where IC± 
are infinite arcs contained in C'±. 

• Case i b) If < rt = TO, the point P is an stable 
node and T does not exists. There exists an hete- 
roclinic sequence of type C_ — > O ^ P or of type 
IC+ — >0 — >P. 



• Case i c) If < 71 < TO < the point P is 
an stable node and T does not exists. There exists 
an heteroclinic sequence of type C_ — > O — > P or 
of type — > O — > P, where k+ is a finite arc 
contained in (7+. 

• Case ii) For ti > and y^n"^ + 1/2 < to < 
^v? + the point T is a stable node and the 
point P is a saddle. For these conditions the hete- 
roclinic sequence is of type C± — > O — > T — > P, 
or of type k+ — > O — > T — > P. 

• Case iii) For ri > and a/71^"+477 < to < \/l + n^, 
the point T is a spiral node and the point P is a sad- 
dle. For these conditions the heteroclinic sequence 
is the same as in the former case. 

• Case iv) For ti > and to > vTTt? the point T 
is a spiral node whereas the point P does not exist. 
The heteroclinic sequence in this case is C_ — s- 
O — >T. 



IV. NORMAL EXPANSION UP TO 
ARBITRARY ORDER 

Normal Form (NF) calculations essentially removes the 
quadratic, cubic, etc., terms that are effectively neglectful 
up to the desired order in the Taylor expansion around 
a non-hyperbolic fixed point, e.g., C±. The terms that 
cannot be removed by using NF calculations, and they 
are the essential degrees of nonlinearity (see appendix 
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TABLE I: Location, existence and deceleration factor of the critical points for m > 0, n > and y > 0. We use the notation 
S — m? — -n? (from reference [3)- 



Name 


X4, 




y 


Existence 


g 


O 











All m and n 


1 
2 


c± 









All m and n 


2 


p 


m 


~n 




5 < 1 


-1 + 3(5 


T 


m 
25 


n 

~2S 


1 


5 > 1/2 


1 
2 



IIV|) . In this section we obtain the normal form expansion 
for nonhyperbolic points in the curve C_ . 

The normal expansion up to order N is given by the 

Proposition IV. 1 Let be the vector field X given by 
SKMIS^) which is C°° in a neighborhood of x* = 
(xj, , y*)"^ G C-. Let m > n > 0, and cc* G R, such 

that A3" = 1 — nx'^ + m^l + x*^ is not integer, then, 

there exist constants Or, br, c^, r > 2 (non necessarily 
different from zero ) and a transformation of coordinates 
X — > y, such that U(M1S\) has normal form 

N 

y'l = J2^ryl + 0{\yf+'), (15) 
y2 = y2 (^1 + ^'■2^1"') + o(|yr+'), (i6) 

2/3 = y3 (^X3 + f2cry[-'^+0{\y\''+'), (17) 

defined in neighborhood of y — (0,0,0). 

Proof. By applying linear coordinate transformations 
one is able to reduce the linear part of X to the form 

/ \ / XI \ 

Xi(x) == 1 X2 = Jx. 

\0 X3 J \X3 J 

By the hypothesis m > n > we guaranteed A3" > 1 for 
all cc* S M. Since, the eigenvalues of J are different and 
J is diagonal; then, the corresponding eigenvectors 

B = |x™ia;^=^x"^e,|TOj G N,^TOj = r,i,j = 1,2,3| 

form a basis of . Since 

Ljx-^e, = {(m-A)-A,}x'"e„ 



the eigenvectors in B for which Am.i = (m • A) — Ai ^ 
form a basis of B^ = Lj(iJ'"). The eigenvectors asso- 
ciated to the resonant eigenvalues, i.e., those such that 
Am,i = 0, form a basis for the complementary subspace, 
G^'of B'' in H''. 

Since Ai = 0, the resonant equations of order r (with 
nil + "I2 + rn^ = r) has unique solution 

^2 + A3m3 = =^ mi = r, m2 — ~ 0, (18) 
m2 + A3TO3 = 1 =^ mi — r — 1, m2 ~ 1, m3 = 0, (19) 
+ A3m3 = A3 mi — r — 1, m2 = 0, ms = 1. (20) 

Equation (j20p has another solution given by mi = r — 
•^3,"T-2 = A3,m3 = 0, provided A3 is integer satisfying 
1 < A3 < r. This solution is discarded by our hypothesis 
on A3. 

Then, \^x\ei,x^^^X2e2,x^^^X3e3} , form a basis for G"' 

in 

By applying theorem lA. 11 we have that there exists a 
coordinate transformation x ^ y, such that (|10m2[) has 
normal form (jlSlllTp where , br and Cr are some real 
constants. I 

The values of of all these constants can be uniquely de- 
termined, up to the desired order, by inductive construc- 
tion in r. The idea is to introduce polynomial coordinate 
transformations x ^ x + hr(x) such that Ljhr(x) = 
Xr(x)-|- higher order terms. The higher order terms of the 
expansion are modified by successive coordinates trans- 
formations. Suppose we apply the r-th transformation, 
for a fixed r. Then, all non-resonant terms of order r, 
are eliminated. The r-th transformation modifies the ex- 
pansion terms of order higher than r, but resonant terms 
of order less than r are not affected. One then removes 
the non-resonant terms of order r -I- 1 by introducing the 
(r-l-l)-th transformation, and so on. In each step one ob- 
tains the desired coefficients. Finally, we can truncated 
the expansion up to order N. 

We see that (|15m7p satisfy the conditions of the exis- 
tence and uniqueness theorem of differential equations. 
Then, there exist a unique solution of (jlSlllTp passing 
through (j/10, 2/20, 2/30 ), at t = 0. By neglecting the error 
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terms one is able to integrate the resulting approximated 
system 



N 



y'l = ^arvl, 

r=2 

y'2 = y2 l^i + ^^r-yr'^ > 
y's = ys ^'^3 + ^CrZ/r^^ 



(21) 
(22) 
(23) 



with initial condition (yi(io), y2(io), 2/3(^0)) = 
(2/10,2/20,2/30)- 

The general solution of (|2m23p is as follows. 

If 2/10 = 0, then yi{t) = 0, 1/2 (0 2/2oe^"^° and y^it) = 
ysoe^^'-''^'"'^ for all t e R. Then, the orbit approach the 
origin as r — > — 00 provided A3 > 0. 

If 2/10 7^ 0, then (|21j|23p can be integrated in quadra- 
tures as 



T -To 



r (fl^rc] dC, (24) 

■^yio \r=2 / 

2/2(0 = 2/2oe"-"° n I 

r=2 L "'^0 

W r 

2/3(<) = 2/306^^^"-"°) n^^P / 



N 



(25) 
(26) 



The yi-component of the orbit passing through 
(2/10,2/20,2/30) at T = To with yio 7^ is obtained by in- 
verting the quadrature ([M)) . 

The other components are given by 



2/2 = 2/20 exp 



2/3 = 2/30 exp 



VI I 
Via 



El2^-C 



r-1 



yi 



J2r=2 O'-C'' 
I + E.I2C.C-' 



-dC 



J2r=2 O-rC 



-dC 



(27) 
(28) 



This orbit does not admit, in general, a prolongation to 
all the time values. If the maximal interval of definition, 
(a, (3) of the solution yi is such that a is finite, then the 
orbits diverges in the yi-direction as t ^ a+. 



A. Treatable case: normal expansion to third order 
for C_ 



In this section we show normal form expansions for 
the vector field (|10m2p defined in a vicinity of C_ ex- 
pressed in the form of proposition IIV.2I The proof is 
given in three steps. First, by applying linear coordinate 
transformations (translation to the origin and similarity 



transformation to reduce the Jacobian matrix to its real 
Jordan form) we reduce the linear part in its simpler 
form. Second, we perform a quadratic coordinate trans- 
formation in order to reduce non-resonant second order 
terms. Third, we perform a cubic coordinate transfor- 
mation that allows to eliminate non-resonant terms of 
third order. In principle it is possible to eliminate all the 
non-resonant terms of all orders. However, we find the 
system computationally treatable to third order. 

Proposition IV. 2 Let be the vector field X given by 
iKMlS^) which is C°° in a neighborhood of x* ~ 



{xZ^,x*,y*) £ C-. Let m > n > 0, and x* G 



such 



that A3 = 1 — nx 



+ my^]~-i-~E^ is not integer, then, 

there exist a transformation to new coordinates x — > z, 
such that diOH jj?)) . defined in a vicinity ofx*, has normal 
form 

z[ ^ 0{\zf), (29) 
2^ - 22 + 0(kh, (30) 

4 = {X3 +C2Zi+C3zf)z3+Oi\zf), (31) 



where C2 — —n- 



and C3 — 



2{l+x 



Proof. 

By using linear coordinates transformations the system 
(fT0llT2)) is reduced to 



x' = Jx + X2(x) +X3(x) 



(32) 



where x stands for the phase vector x = {xi, X2, X3) , J 
stands for the Jordan Form of the matrix of derivatives 



/O 
1 



1- nx*^ + mJl + xy 



(33) 



X2(X) 



-'^(l,l,0),l2^1 ^2 + -''^(0,0,2) 42:3 
-^^(2, 0,0), 22^1 + -^(0,2,0), 22^2 + -'^(0, 0,2) ,22^3 

\ -'^(i,o,i),3a;ia;3 +^(0,1,1), 3a;2a;3 / 



(34) 

where the coefficients Xm,j: are displayed in table [Til for 
the allowed m. 

The third order terms become: 



X3(x) 



/ -'^(3,o,o),ia;i + ^(1,2,0), 12^12^2 + -'^(i,o,2),ia;ia;§ ^ 
-'^(2,i,o),2a;ia;2 + -^^(0,3,0), 22^2 + -'^(o,i,2),2a;2a;i 



V -'^(2.o,i).3a;ia;3 + ^(o.2,i).32;la;3 + ^(o,Q.3),3a;i / 

(35) 

with the coefficients defined as in table lllll 

Second step: simplifying the quadratic part 
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TABLE II: Coefficients of the vector field X2(x). 



m 


Xin,l 


Xra,2 


Xtn,3 


(2,0,0) 





2(.;2+i) 





(1,1,0) 










(1,0,1) 










(0,2,0) 





3 





(0, 1, 1) 










(0,0,2) 


mx* -^1 + 2 -n(^l + X* ^) 


ix; (-2771^1 + a;;2 + 2nx*^ - l) 






TABLE III: Coefficients of the vector field X3(x). 



m 




Xm,2 


^m,3 


(3,0,0) 


1 

2(x*2+l) 








(2,1,0) 





1 

2(x;2+i) 





(2,0, 1) 








1 

2(x;2+i) 


(1,2,0) 


1 

2xJ^ 








(1,0,2) 


1 
2 








(0,3,0) 





1 

2^ 





(0,2, 1) 








1 

2^ 


(0,1,2) 





1 

2 





(0,0,3) 








1 

2 



TABLE IV: Eigenvalues of L(2)j -.H^^H^. 



m Am,l Am,2 Am,3 



(2,0,0) 


-1 




(1,1,0) 1 






(1,0,1) 







(0,2,0) 


1 




(0, 1, 1) 




1 


(0,0,2) 2 (1 - ns; + my'l + a;;2) 


1 - 2nx*^ + 2my/l + x*^^ 





By the hypotheses the eigenvalues Ai = 0, A2 = its eigenvectors form a basis of . The linear operator 
1, A3 = 1 — nx* ^ + m-y/l + a;* ^ of J are different. Hence, j^Cz) . ^2 _^ jj2 



has eigenvectors x^e^ with eigenvalues A^^ i = toiAi + 
m2A2 + A3m3 — Aj, i = 1,2,3, mi, 7712,77x3 > 0, mi + m2 + 
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1713 = 2. The eigenvalues Am,i for the available m,i are 
displayed in table IIVI To obtain the normal form of (|29f 
I3ip wee must look for resonant terms, i.e., those terms 
of the form x^e^ with m and i such that Am.i = for 
the available m,i. Only one term is resonant of second 
order: A(i,o,i),3 = ^ 222/12/363. 
The required function 



to eliminate the non-resonant quadratic terms is given by 

/ #^^yiy2 + #^^2/| \ 



h2(y) 



^(2, 0,0), 2 2 _(_ ^(0,2,0), 2 2 _|_ ^(0,0, 2), 2 2 
Aro n 01 9 "1 A/n 9 01 9 "2 A/n n 91 9 ""^ 



V 



^(0,2,0), 2 



A, 



(0,1,1), 3 



2/22/3 



The quadratic transformation 

X ^ y + h2(y) 



(36) 
(37) 



with h2 defined as in (|36p is the coordinate transforma- 
tion required in theorem lA.il By applying this theorem 
we prove the existence of the required constant C2. 

To finish the proof, let us calculate the value of C2 . 

By applying the transformation ([57)1 the vector field 
([5^ transforms to 

y' = Jy - Lf h2(y) + X2(y) + X3(y) + 0(|y|^), (38) 
Since 

- L^^^h2(y) + X2(y) = ^(i,o,i),3yi2/3e3, (39) 
we have 

y' = Jy + X(i,o,i),32/i2/3e3 + X3(y) + 0{\y\^) (40) 

, i.e., C2 = A:(i,o4),3 = -n+ ^/^rf- 

The vector field X3(y) introduced above has the coef- 
ficients: 



A'(2,0 



(2, 0,1), 3 



1 2(a:;2, 



X 



(1,2:0), 1 - —2' 



-'^(1,0,2),1 — ^ 
-'^(1,0,2),2 = 



2(^2 + TO ((5 + m _ 1)) _ na;* [2m(5 + + l 



(A3 



„ ^ „ fa;* -2f2((5^-l)ri+ri))m 

2a;* + ^ — 5 + " 



(2na;* 



2A; 



(0,3,0), 2 



X, 



(0,2,1), 3 



-^^(0,1,: 



(-m(5 + nx* - 2) (-2m(5 + 2na;* - 3) 
5 (4(52a;* TO^ + U/^irri^ + A2TO -|- nJAa) 



(0,1, 2), 2 



2 ((^2 - 1) + x;[4m^ + 3]n - m(5(2m(5 + 3) - 3, 

[na;* - 2] - 5, 



X, 



(0,0, 3), 3 



-4nx^ [na;^ 



r 



(41) 



where S = \/^" 
, 2 



+ 1, Ai = 2x; - n (3x;2 + 1) , 

A2 - 4n(-4x;2 + n(3a;;2^2)x;-2) + 5x;, A3 = 
-4na;* [na;* - 2] - 5. 



Third step: simplifying the cubic part 



After the last two steps, the equation is trans- 
formed to 

y' = Jy + 022/12/363 + X3(y) -t- 0(|y|4). (42) 

There is only one term of order three which is resonant 
(see table El): A(2,o,i),3 = ^ 032^2383. 

As in the last step, in order to eliminate non-resonant 
terms of third order we will consider the coordinate trans- 
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TABLE V: Eigenvalues of L^^^ : ^ . 

m Am.l Am, 2 Am, 3 

(2,0,1) - - 

(1,2,0) 2 

(1,0,2) 2 (m^l + a;;2 _ + 1) 2m^l + a;;^ _ 2nx*^ + 1 

(0,3,0) - 2 

(0,2,1) - - 2 

(0,1,2) 2m^l + a;;2 -2nx*^ + ?, 2 {m,/TTx^ - nx*^ + l) 

(0,0,3) - - 2(myrT^-na;; + l) 



formation y — > z given by 

y = z + h3(z) 

where 



(43) 



is defined by 



h3(z) 



where Ar 
erator 



/ ^(l,2,0),l ^ ^2 

A(i,2,o),i 2 



A 



(1,0, 2), 2 
\ A(o,2,l),3 



HI. 0.2), 
'^(0.3.0). 2 3 

^2" 



(0,1,2),1 ^2 \ 



X. 



(0,2,1), 3 ^2^ I ^ 

2:3 2:3+ 



^(0,3,0), 2 
(0.0. 3), 3 „3 



(0.1. 2). 2 ^ 2 
72:2^3 



(0,0. 3), 3 



(44) 

are the eigenvalues of the operator linear op- 



associated to the eigenvectors x^e^. In t able fVl are shown 
these eigenvalues. The associated eigenvectors form a 
basis of (the space of vector fields with polynomial 
components of third degree) because the eigenvalues of J 
are different. 

The transformation (|43p is the required by theorem 
(|A.1[) . By using this theorem we prove the existence of 
the required constant C3 . To find its we must to calculate 

-Lf h3(z)+X3(z), 

which it is equal to 

where e^, i — 1,2,3, is the canonical basis in M". Then, 
m n-x* 

C3 = 



Observe that the transformation h3 does not affect the 
value of the coefficient of the resonant term of order r = 
2. Then, the result of the proposition follows. ■ 



1. Unstable manifold to third order 

For A3" > 0, the origin has a 2-dimensional local un- 
stable manifold tangent to the plane 2:2-2:3 at given by 

WilM = {(21, 22, 23) e : 21 = h{z2, 22), 

Dh(0) = 0, (22,23)^ <<5}, (45) 



where /i : — > R is a C function and 5 > is small 
enough. 

Using the invariance of W^;"^(0) under the dynamics of 
(|29ll3ip we obtain a quasi-linear partial differential equa- 
tion that h must satisfy: 



AA(/l(22,23)) = 0(|z|4) (46) 

where we have defined the differential operator 



dh dh 
N{h{z2,z:i)) = 22 — +23 (c3/i^ +C2/1 + A3-) — . (47) 
0Z2 0Z3 

This system is solved up to order 0(|z|^). 
Assuming that Wf^^{0) is C* (M^) 

we can expressed it 

as the graph 



h (22, 23) = /i3,o2| + ft-2,o2| -I- 23/12,12! -I- 23/11,122 -I- 

23/11,222 + 2|/io,2 + 2f/io,3 + 0(|z|4).(48) 

This function satisfy the tangentiality conditions 
h{0, 0) = d,,h{0, 0) = d,,h{0, 0) = (0, 0). 

By substitution in the differential equation 
J\f {h (22, 23)) = and by discarding the error terms we 
find 
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3/i3,o4 + 2/12,0^2 + '^2423 (A;^ + 2) + /ii,iZ3 (AJ + l)z2 + hiazl [2\^ + l]z2 + 2z^ho^2>^3 + Szf/io.sAg = (49) 



r 



Using the condition A3 > we get h^^ = /i2,i — 
^2,0 = ^1.2 = ^1,1 = ^0,3 = ^0,2 ~ 0. Then, the un- 
stable manifold of the origin is, up to order 0(|z|'*), 
W^^ocW = {(^1,^2,^3) e : z, = 0,z| + z| < ,52} 
where (5 is a real value small enough. Therefore, the 
dynamics of (|29ll3ip restricted to the unstable mani- 
fold, is given, up to order Odzj*), by zi = 0, Z2{t) = 

e'^Z2o, Z2{t) = e^^ ^2^30, where Z2q + z'^q < S^. This means 
that linv„,_oo {zi{t), Z2{t), Z3(r)) = (0,0,0) . Then, the 
origin is the past attractor for an open set of orbits of (|29f 
I5T|) . In the original coordinates this means that the past 
asymptotic dynamics of an open set of orbits of the sys- 
tem can be accurately approximated by a massless quin- 
tom model located at C_. Let us investigate the center 
manifold. 



2. Center manifold to third order 

For A;^ > 0, the origin has a 1-dimensional local un- 
stable manifold tangent to the plane zi at given by 



WU0)={{zi,Z2,Z3)eM.' 

Df{0) = 



Z2 = f{zi), Z-i = 5(zi), 
0,i?5(0) = 0,|zi| <,5}(50) 

M are C' functions of z\ 



where / : M ^ M and g : M - 
and (5 > is small enough. 

Using the invariance de under the dynamics 

of (|29ll3ip we obtain a system of two ordinary differential 
equations that / and g must satisfy: 

/(zi)-/'(zi)0(|zn=0(|z|4),(51) 
(A3- + C2Z1 + C3Z?) g(zi) - a\z^)0 (|z|4) = O (|z|'i) (52) 

In order to obtain an expression for the center mani- 
fold of the origin up to order O (|z|'*) we need to solve 
equations (|5m52p up to the same order. Assuming that 
^ioc(^) expressed as 



The last implication follows from the fact that A3 > 0. 
Then, the center manifold of the origin is, up to order 
0(|z|4), W^f„,(0) = {(zi, Z2,z^) e R3 : Z2 = Z3 = 0, |zi| < 
(5} where (5 is a real value small enough. 

Let us describe the center manifold, up to the pre- 
scribed order, as a graph in the original variables. 

We have 



2^x;2 + 1 

(57) 

+ zi+0(|zi|4), (58) 



X2 — -^ip 

X3 = y = {\zi\ 



2x^2 _^ 2 



(59) 



By taking the inverse, up to fourth order, of ([58]) we 
have the expression for zi: 



z\ 



* 2 3 



Substituting (|5D|) in (|57ll59p we have that the center man- 
ifold of the origin is given by the graph: {(xi,X2,X3) G 



2(^.2+1)5/2 

= 



X2 



2(^-2 + ij3/2 ^ 

X2I < (5} for (5 > small 



0(|X2|4),X3 

enough. 

In conclusion, the analysis up to fourth order does not 
allows to classify the global past asymptotic dynamics of 
quintom model. Since the numerical experiments in Q 
suggest that there is an open set of orbits that tends to 
infinity, then it is worthy to investigate the dynamics at 
infinity. 



V. ANALYSIS AT INFINITY 



/(zi) =az2 + 5^3^0(|zi|4), 
g(zi) =cz2 + dz3 + 0(|zi|4). 



(53) 
(54) 



The representation (|53ll54p satisfy the tangentiality con- 
ditions /(O) = .g(0) = /'(O) = g'(0). By substitution of 
(|53ll54p in (|5m52p and neglecting error terms we get 



az\ 



0(|zi|4) a^6==0, 
(cc2 + d\l) z\ + cX^zl = O (Izil-*) 



cAg = CC2 

c = d = 0. 



(55) 



(56) 



In this section we investigate the dynamics at infinity 
by using the central Poincare projection method (sij . 

To obtain the critical points at infinity we intro- 
duce spherical coordinates (p is the inverse of r = 

I, -f x^ -I- y^, then, p ^ as r — > 00): 



Xfh = — sin 9i cos 6*2 , 
P 



y = — smWi sin6'2, 
P 



Cm = — cose/i 
P 



(61) 
(62) 
(63) 
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TABLE VI: Location and existence conditions for the critical points at infinity. 



Name 


01 


62 


Existence 







2 


always 


A* 


TV 


2 


always 


-'3 


IT 

4 


±cos-i(-f) 


-TT < ±cos"^ (-^) < 7r,n / 


p± 


3ir 
4 




-TT < icos^^ (^) < 7r,n / 


P5 







< < TT 








< < TT 



TABLE VIL Stability of the critical points at infinity. We use the notation 5 — m? — r? and A* = n cos 61 ± msin6'J'. 

Name (Ai,A2) p' Stability 

P± {-n,n) 
Pt (-n,n) 

p± / V2S S \ 

^3 \ n ' yf2n) 



p± ( _^/M. 

P^ (0,A+) 
fe (0,A-) 



where < 6*1 < tt and — tt < 6*2 < tt, and < p < cxd. (|10m2p . can be written as 

Defining the time derivative /' = yOd//dr, the system 

I 

p' = 1 (cos^ 6*1 - cos(26i2) sin^ 6I1) + 2n cos 6*1 sin^ 6I1 sin^ B2P + O (p^) . (64) 

and 

Q'^=n cos(26'i) sin 6*1 sin^ 6*2 - cos 6*1 sin 6*1 sin^ 62P + O [p^) , 

62 — (n cos 01 cos 02 + TOsin^i) sin 02 — cos 02 sin 02/9 + O (p^) . (65) 

I 

Since equation (|64p does not depends of the radial com- the compatibility conditions 
ponent at the limit p — s- 0, we can obtain the critical 
points at infinity by solving equations (j65p in the limit 
p 0. Thus, the critical points at infinite must satisfy 



> 

> 

> 0, S <0 
<0, 5>0 

> 0, S <0 
< 0, S>0 



< 0, f < 



4 



> 0, otherwise 



<0, 



< 



4 - '^1 < 4 
> 0, otherwise 



3ir 



saddle 
saddle 

source if n < 0, n < m < — n 
saddle otherwise 

source if n > 0, ~n < m < n 
saddle otherwise 

nonhyperbolic 

nonhyperbolic 
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cos(26'i)sin6'i sin^ 6*2 = 0, 
(n cos cos 6*2 + TO sin 9i ) sin 02 = 0. (66) 

First, we examine the stability of the pairs {6\, 62) sat- 
isfying the compatibihty conditions (|66p in the plane 61- 
02, and then, we examine the global stability by sub- 
stituting in ((64|) and analyzing the sign of p'(02j^2)- 
table IVll it is offered information about the location and 
existence conditions of these critical points. In table IVlJ 
we summarize the stability properties of these critical 
points. 

The cosmological solutions associated to the critical 
points and have the evolution rates 4''^/V = 
0, = and H/ip = p/V6 ^ 0. ^ These solutions 
are always saddle points at infinity. The critical points 
P^ and P^ are sources provided n < 0,n < m < —n or 
n > 0, — n < m < n, respectively. They are saddle points 
otherwise. The associated cosmological solutions to P^ 

have the evolution rates ijy^/V — , 4>/^ = —m/n, 

and H/(j) = —np/{\/3m) — > 0, and H/ip = p/VS 0, 
whereas the associated cosmological solutions to P^ have 

the evolution rates 4>'^ /V = rfi^-m^ ■> i'l'f = —m/n, and 
H/^ = np|{^/lm) 0, and H/ip = -p/Vi 0. 
The curves of critical points P5 and Pg ^ire nonhyper- 
bolic. The associated cosmological solutions have ex- 
pansion rates (vahd for 9^ ^ 7r/4) V/(p'^ — 0,(p/(p = 
t&ne\,Hl^p = psece\/^/& ^ 0, and V/4>'^ = 0,<j)/ip = 
— tan 9*, H/ip — p sec 9\ /\/6 — > 0, respectively. 

Concluding, although there exists unbounded orbits 
towards the past, by examining the orbits at infin- 
ity, we get that the sources satisfy the evolution rates 
(jy^ /V ~ , 0/<y9 ~ —m/n, with H/cf) — > 0, provided 

n < 0, n < m < — n or n > 0, —n < m < n. 

VI. CONCLUDING REMARKS 

In this paper we have investigated the past asymptotic 
dynamics of quintom cosmologies (which is closely related 
with the behavior of typical orbits near an hyperbolae of 
critical points) by using the Normal Form Theorem as 
a tool. We have developed the general (arbitrary order) 
normal expansion for the quintom model. This expan- 
sion is the minimal in the sense that the terms that are 
involved are the "essential" degrees of nonlinearity. The 
other higher order terms are removed by making properly 
coordinate transformation. By integrating the resulting 
system one is able to construct unstable and center man- 
ifolds. We have explored the structure of unstable and 
center manifold of the origin up to fourth order. From the 



structure of the unstable manifold we see that that the 
past asymptotic behavior of quintom cosmologies, with 
exponential potentials, is given by a massless scalar field 
cosmology for an open set of orbits. We provide here an 
approximated formula for the center manifold. By exam- 
ining the center manifold we find that, up to fourth order, 
the curve of critical points C_ is not the past asymptotic 
global attractor. 

To get more accuracy we need to deal with the expan- 
sion up to arbitrary order iV > 3. However, it is very 
difhcult to get the expansion coefficients for N > 3. The- 
oretically, if the solution of (pij) for yi admits a prolonga- 
tion to —00 < r < 00, (for instance, if yio = 0), then the 
trajectory passing through (yio, y20, 2/30) at to approaches 
the origin as r ^ —00, provided A3 > 0. Otherwise, if 
the maximal interval, (a,/?), has a finite, then, the or- 
bits diverges in a finite time has r — > a^. This fact is 
supported by the numerical simulations in reference 
(it appears that an open set of orbits escape to infinity 
into the past, which means that the local sources could 
not be the past attractor). 

By examining the dynamics at infinity, using the cen- 
tral Poincare projection method, we get that the sources 
satisfy the evolution rates ip^ /V ~ , (j>/(p ^ —m/n, 

with H an infinitesimal of 0, provided n < 0,n < to < 
— n or n > 0,— n < m < n. These results complete the 
past asymptotic analysis of quintom cosmologies. 
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APPENDIX A: NORMAL FORMS FOR VECTOR 
FIELDS 

In this section we offer the main techniques for the 
construction of normal forms for vector fields in K". We 
follow the approach in [l^l • 

Let X : M" M" be a smooth vector field satisfying 
X(0) = 0. We can formally construct the Taylor expan- 
sion of X about 0, namely, X = Xi + X2 -I- ... -I- X^ -|- 
0(|x|'"'+^), where X^. G H'^, the real vector space of vector 
fields whose components are homogeneous polynomials of 
degree r. For r = 1 to fc we write 

Xr(x) — J2mi = l ■ ■ ■J2m,^ = l ^j=l ^m,jX™ej, 

E^m^=r, (Al) 

Observe that Xi = DX(0)x = Ax, i.e., the matrix of 
derivatives. 

The aim of the normal form calculation is to construct 
a sequence of transformations which successively remove 
the non- linear term X^, starting from r — 2. 

The transformation themselves are of the form 



^ Do not confuse p with the matter energy density, the latter de- 
noted by pm- 



x = y + h,r(y). 



(A2) 
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where G , r > 2. 

The effect of (|X2|) in Xi is as fohows [13]: Observe that 
X = 0(|y|). Then, the inverse of (jA2| takes the form 

y = x-h,,(x) + 0(|xr+i). (A3) 

By applying total derivatives in both sides, and assum- 
ing x' = Ax + Xr (x) , we find 

y' = Ay - LAh,(y) + X,(y) + 0{\y\^+') (A4) 

where La is the linear operator that assigns to h(y) £ 
the Lie bracket of the vector fields Ay and h(y): 

La : i?'' ^ H'' 

h ^ LAh(y) = Dh(y)Ay - Ah(y). (A5) 

Both La and X^ G H^, so that the deviation of the 
right-hand side of (jA4p from Ay has no terms of order 
less than r in |y|. This means that if X is such that 
X2 = ...Xr-i — 0, they will remain zero under the 
transformation (jA2p . This makes clear how we may be 
able to remove X^ from a suitable choice of X^. 

The proposition 2.3.2 in [l^] states that if the inverse 
of La exists, the differential equation 

x' Ax + X,,(x) + 0(|x|''+i) (A6) 

with Xr G H'', it is transformed to 

y'-Ay + 0(|yr+i) (A7) 

by the transformation (jA2p where 

h.(y) =LA-^X,(y) (A8) 

The equation 

LAh,(y) = X,(y) (A9) 

is named the homological equation. 

If A has distinct eigenvalues A.;, i — 1,2,3, its eigen- 
vectors form a basis of M". Relative to this eigenbasis, 
A is diagonal. It can be proved (see proof in [27l|) that 
La has eigenvalues Am,i = m • A — = J2j ^T^j^j ~ 
with associated eigenvectors x'^e^. The operator, La~^, 
exists if and only if the Am,i 7^ 0, for every allowed m 
and i = 1 . . .r. 

If we were able to remove all the nonlinear terms in this 
way, then the vector field can be reduced to its linear part 

x' - X(x) ^ y' = Ay. 

Unfortunately, not all the higher order terms vanishes by 
applying this transformations. Particularly if resonance 
occurs. 



The n-tuple of eigenvalues A = (Ai, . . . , An)"^ is reso- 
nant of order r (see definition 2.3.1 in ^]) if there exist 
some m — (mi, m2, . . . mn)^ (a n-tuple of non-negative 
integers) with mi -I- m2 -|- . . • m^ = r and some i = 1 . . . n 
such that Ai = m • A, i.e., if Am.i — for some m and 
some i. 

If there is no resonant eigenvalues, and provided they 
are different, we can use the eigenvectors of A as a basis 
for . Then, we can write as 

hr.(x) = ^ h^^i-X^Bi 

m,'i,^ rrij—r 

and any vector field X G H'^ as 
X(x) = 

m,z, — ^ 

where m = (mi, m2, . . . , mn)^ , x"* = a:™^ x™^ . . . a:™" 
and Si, i — 1, . . .n stands for the canonical basis in R". 
If the eigenvalues of A are not resonant of order r, then 

This gives explicitly in terms of X^. 

In the case of resonance occurs, we proceed as follows. 
If A can diagonalized, then the eigenvectors of La form 
a basis of . The subset of eigenvectors of La with non- 
zero eigenvalues then form a basis of the image, B^, of 

under La. It follows that the component of X^ in 
can be expanded in terms of these eigenvectors and 
chosen such that 

to ensure the removal of these terms. The component, 
Wr, of Xr lying in the complementary subspace, G"", of 
B^ in will be unchanged by the transformations x = 
y -I- h.r{y) obtained from B^ . 
Since 

X,(y+h,+fc(y)) = X,(y)+0(|y|'-+^-+i), r > 2, fc = 1, 2, . . . , 

these terms are not changed by subsequent transforma- 
tions to remove non-resonant terms of higher order. 
The above facts are expressed in 

Theorem A.l (theorem 2.3.1 in [27]) Given a 
smooth vector field X(x) on R" with X(0) = 0, there 
is a polynomial transformation to new coordinates, y, 
such that the differential equation x' = X(x) takes the 

form y' = Jy + I]^i"Wr(y) + 0(|y|^+^), where J is 
the real Jordan form of A. = DX(0) and w,, £ , a 
complementary subspace of H'^ on B^ = LA(-ff'^). 
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